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ABSTRACT 


Two  objective  studies  are  made  of  geostrophic  relative  vorticity 
computed  from  the  height  field  and  of  observed  relative  vorticity  com¬ 
puted  from  wind  components  at  500  mfe. 

In  the  first  study  an  assumption  was  made  of  linear  wind  variation 
such  that  observed  vorticity  on  a  600 km  scale  was  assumed  to  be  re¬ 
presentative  of  the  300  km  scale  on  which  the  geostrophic  vorticity  was 
computed.  The  very  poor  agreement  of  the  vorticities  computed  by  this 
method  indicated  that  this  assumption  was  erroneous. 

A  second  study  was  made,  therefore,  of  600  km  scale  observed 
vorticity  and  of  300km  and  600  km  scale  geostrophic  vorticity.  The 
comparison  of  geostrophic  vorticity  of  different  scales  indicated  the 
importance  of  grid  size  in  the  computations.  The  comparison  of  ob¬ 
served  and  geostrophic  vorticity  of  the  same  scale  indicated  small 
though  significant  differences. 
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A  COMPARISON  OF  RELATIVE  VORTICITIES 
COMPUTED  FROM  GEOSTROPHIC  AND  OBSERVED  WINDS 

I.  INTRODUCTION 

Within  recent  years  the  vorticity  concept  has  come  into  increasingly 
greater  use  in  both  synoptic  and  dynamic  meteorology.  In  most  studies 
of  vorticity  the  geostrophic  approximation  has  been  used  to  compute  the 
relative  vorticity  from  the  height  field  of  a  constant  pressure  chart  in 
preference  to  the  use  of  the  upper-air  wind  field.  Many  theoretical  and 
practical  reasons  have  been  proposed  both  to  justify  and  to  qualify  the  use 
of  this  technique.  In  numerical  weather  prediction,  for  example,  the  use 
of  the  geostrophic  assumption  is  advantageous,  not  only  as  a  filter  to  re¬ 
move  certain  undesirable  wave  solutions,  but  aB  a  convenient  means  of 
dealing  directly  with  functions  of  a  scalar  height  field  rather  than  re¬ 
ducing  wind  fields  into  the  scalar  magnitudes  of  components  or  making 
separate  computations  of  wind  shear  and  curvature.  Recent  studies  of 
cyclogenesis  (Spar,  et  al,  ,  1955)  have  been  based  upon  the  initial  deter, 
mination  of  geostrophic  vorticity  before  the  computation  of  the  advective 
and  generative  terms  of  the  vorticity  equation  at  various  upper  levels 
above  regions  of  cyclogenesis.  The  use  of  the  vorticity  concept  in  meteor¬ 
ology  has  become  so  increasingly  important,  and  the  use  of  the  geostro¬ 
phic  assumption  is  so  widespread,  that  an  estimate  of  the  validity  of  the 
geostrophic  technique  is  of  considerable  concern. 
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Much,  study  has  already  been  devoted  to  the  deviations  of  the  geo- 
strophic  from  the  observed  wind,  and  more  especially  in  recent  years 
to  the  vorticity  deviations  which  naturally  arise  from  such  differences. 
Charney  (1948)  laid  much  of  the  basis  of  current  computational  techniques 
in  numerical  weather  prediction  by  reasoning  on  theoretical  grounds  that 
the  geostrophic  deviation  is  negligible  for  those  primary  large-scale  per¬ 
turbations  where  the  mean  horizontal  distance  between  trough  and  wedge 
is  of  the  order  of  1000  km.  His  conclusion  was  based  on  the  assumption 
of  typical  values  for  the  orders  of  magnitude  of  the  mean  horizontal  dis¬ 
tance  between  trough  and  wedge,  the  mean  vertical  distance  between  points 
at  which  velocity  components  assume  extreme  values,  the  mean  magnitudes 
of  the  horizontal  and  vertical  velocity  components,  and  the  mean  speed  of 
propagation  of  the  system.  Charney  was  thus  able  to  show  that  the  hori¬ 
zontal  acceleration  of  a  1000  km  perturbation  is  one  order  of  magnitude 
less  than  that  of  the  horizontal  Coriolis  force. 

Among  the  early  synoptic  investigations  of  geostrophic  wind  devia¬ 
tions  was  that  by  Houghton  and  Austin  (1946),  who  studied  the  deviations 
of  the  geostrophic  wind  at  10,  000  feet.  Thiyi  measured  the  geostrophic 
wind  components  from  the  spacing  between  isobars  and  obtained  the  actual 
wind  components  from  PIBAL  and  RAWIN  observations.  By  then  subtract¬ 
ing  the  geostrophic  from  the  observed  wind  component,  they  obtained  the 
geostrophic  deviations,  u'  and  v'.  For  619  cases  they  computed  a  mean 
absolute  value  of  u1  and  v'  as  7.0  mph,  compared  to  mean  absolute  values 
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of  the  geostrophic  components,  Ug8  =  18.7  mph  and  Vg8  =  12.7  mph. 

The  magnitude  of  the  mean  vector  difference  between  the  actual  and  the 
geostrophic  wind  was  10.  8  mph.  The  deviations  u1  and  v'  were  accepted 
as  real  since  their  isopleths  formed  systematic  patterns. 

Neiburger  et  al.  (1948)  performed  geostrophic  and  gradient  wind 
computations  for  all  wind-reporting  stations  on  two  700-mb  charts  (135 
cases  in  all).  A  variability  of  about  25  percent  was  found  in  the  compu: 
tations  made  by  different  individuals  in  the  same  cases.  The  computed 
geostrophic  wind  speeds  differed  by  about  35  percent  from  the  observed 
wind  speeds.  A  study  of  deviations  as  a  function  of  curvature  produced 
no  significant  correlation,  although  computed  winds  were  predominantly 
greater  than  the  observed  when  the  contour  curvature  was  anticyclonic 
or  neutral,  and  more  evenly  distributed  between  greater  and  smaller 
values  when  the  curvature  was  cyclonic. 

Godson  (1950)  studied  618  cases  at  10,000  feet  for  a  3-day  sequence 
over  the  North  American  continent.  He  computed  a  mean  absolute 
angle  of  14.4°  ±  0.  7°  between  the  reported  and  the  geostrophic  winds. 
The  mean  absolute  deviation  of  the  speeds  was  approximately  7  mph. 

For  latitudes  north  of  30°  N  and  for  a  range  of  geostrophic  wind  speed 
between  20  and  45  mph,  the  difference  between  observed  and  geostrophic 
wind  is  approximately  20  percent  of  the  true  wind  speed. 

Reed  (1951)  considered  the  effect  of  using  geostrophic  winds  for 
vorticity  computations  at  4  and  10  thousand  feet.  He  computed  relative 
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vorticity,  J  ,  from  the  following  equation,  which  expresses  the  shear  and 
curvature  terms  in  natural  coordinates: 

3  =-  8v/Sn  +  v3c/9s  , 

where  v  =  velocity  of  the  wind,  n  =  direction  normal  to  the  streamlines, 
a  =  wind  direction  in  radian  measure,  and  s  =  distance  in  the  streamline 
direction.  Graphical  addition  was  used  to  add  the  separate  shear  and 
curvature  components.  From  the  results  of  his  study,  Reed  concluded 
that  geostrophic  vorticity  may  be  substituted  for  the  actual  whenever  the 
spacing  and  curvature  of  the  contours  can  be  determined  with  a  minimum 
of  ambiguity. 

Newton  (1954)  studied  the  effect  of  contour  analysis  on  the  computed 
vorticity  field.  Geostrophic  vorticity  calculations  were  found  to  be  very 
sensitive  to  details  of  analysis  since  second-order  differentials  of  the 
height  field  are  used.  He  suggested  that  arbitrary  variations  in  the 
vorticity  field  be  minimized  by  smoothing.  Newton  also  pointed  out  that 
the  error  in  computed  geostrophic  wind  varies  inversely  with  the  dis¬ 
tance  between  stations,  for  a  given  height  error. 

Hovmdller  (1952)  considered  the  effect  of  temperature  measurement 
errors  in  the  following  manner:  if  we  have  two  stations,  A  and  B,  the 
relative  error  of  the  geostrophic  wind  perpendicular  to  AB  can  be  shown 

*AVA  JaT A'  aXB>  ZA 
<V„  "  ZA-  %  TA 


to  be 
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where  (ATA-  ATg)  is  the  difference  of  the  mean  temperature  errors 
of  the  station  soundings  for  the  entire  layer  below  the  level  considered. 
At  500mb,  if  we  assume  as  typical  values:  ZA  =  5700m,  TA  =  25Q°A, 
(ZA.  Zg)  =  5&m,  and  (ATA  .  ATB)  =  1°A, 


(Vg)n 


4&o/0 


On  the  basis  of  these  considerations  Hovmidller  declared  that  it  was 
inadvisable  to  use  unsmoothed  height  values  for  computing  geostrophic 
vorticity;  in  fact,  he  stated  that  in  many  cases  at  500  mb  the  observed 
wind  normally  may  be  a  better  approximation  to  the  geostrophic  wind 
than  the  value  derived  from  the  height  field.  He  concluded  by  stating 
a  preference  for  the  following  method  to  compute  relative  vorticity 
from  the  observed  winds 

'S  =  Av/Ax  -  Au/Ay 

where  Au  and  Av  are  finite  differences  of  the  wind  velocity  components 
as  measured  along  intervals  Ax  and  Ay  which  are  centered  at  the 
origin  of  a  Cartesian  coordinate  system. 

Aubert  (1953)  showed  that  the  use  of  the  geostrophic  wind  approxi¬ 
mation  in  relative  vorticity  calculations  at  500  mb  is  equivalent  to 
using  analyzed  wind  fields  having  a  standard  error  of  9  mps  (these 
results  were  based  on  calculations  made  over  5-degree  "squares"). 
Since  this  is  considerably  larger  than  the  value  of  the  standard  wind  . 
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error  obtained  by  Rapp  (1952),  he  concluded  that  the  relative  vorticity 
can  beat  be  obtained  from  observed  winds.  Aubert  analyzed  500 mb  maps 
for  four  synoptic  periods,  analyzed  the  wind  fields  by  the  Bjerknes  isogon- 
isotach  technique,  the  height  contours  by  the  usual  technique  of  subjective 
interpolation,  and  in  each  analysis  made  visual  reference  to  the  other. 

The  geostrophic  wind  at  the  grid  points  was  computed  over  a  length  incre¬ 
ment  of  2.5  degrees  of  latitude.  The  relative  vorticity  was  then  computed 
from  the  wind  values  as  the  average  circulation  per  unit  area  over  5  and 
10  degree  "squares". 

For  the  5  degree  "square",  Aubert  obtained  a  linear  correlation 
coefficient  r  =  0.529  for  242  cases,  a  standard  deviation  of  relative  vorti¬ 
city  of  observed  wind,  or  =  3.81  x  lO'^  sec**,  a  standard  deviation  of 
relative  vorticity  of  geostrophic  wind,  =  4.59  x  10*^  sec**,  and  a  stand¬ 
ard  error  of  estimating  the  relative  vorticity  of  the  observed  wind  from 
that  of  the  geostrophic  wind,  s  =  1.23  x  10*5  sec-*.  For  a  10-degree 
grid,  Aubert  obtained  these  values:  r  =  0.773,  <r  =  2.80  x  10*5  sec*l,  . 

Og  =  2.24  x  10"^  sec"*,  and  s  =  1.78  x  10*5  sec-*.  The  latter  results  show 
the  effect  of  greater  space  smoothing  obtained  by  quadrupling  the  area 
over  which  the  circulation  integral  was  evaluated. 

In  connection  with  a  vorticity  study  of  east  coast  cyclogenesis, 

Spar  et  al.  (1955)  compared  vorticity  values  computed  by  the  wind  com¬ 
ponent  method  with  those  obtained  by  the  geostrophic  method.  500  mb 
data  for  three  upper-air  observation  times  are  used,  and  vorticity 
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values  on  a  300  km  scale  were  computed  from  analyzed  charts  of  height 
contours  and  wind  components.  A  comparison  of  approximately  225 
absolute  vorticity  values  revealed  a  close  agreement  between  the  means, 
only  slight  differences  between  respective  maxima  and  minima,  and  a 
close  similarity  between  analyzed  vorticity  patterns.  The  contour  analysis 
was.  However,  drawn  to  fit  the  wind  data  so  that  the  two  computations  were 
not  based  upon  independent  data. 

A  serious  criticism  of  many  previous  studies  of  geostrophic  vorti¬ 
city  deviations  is  their  basis  upon  synopitcally  dependent  analyses  of  winds 
and  contours,  as  a  result  of  which  both  real  and  unreal  differences  are 
smoothed  out.  The  great  advantage  of  a  smoothing  technique  is  the  cor¬ 
rection  applied  to  erroneous  and  unrepresentative  observations;  the  dis¬ 
advantage,  however,  is  the  large  degree  of  subjectivity  involved  in  the 
many  decisions  of  accuracy  and  representativeness  which  are  necessarily 
made  in  the  analysis  of  a  constant  pressure  chart.  Miller  (1947)  remarked 
that  when  the  number  of  pressure  and  wind  data  is  not  infinite,  it  is 
always  possible  to  construct  isobars  such  that  every  observed  wind  is 
equal  to  the  geostrophic  wind.  Since  such  a  wide  degree  of  discretion 
is  available  to  the  analyst,  the  degree  of  correspondence  between  geostro¬ 
phic  and  observed  winds  and  velocities  is  as  variable  as  the  ability,  ex¬ 
perience  and  synoptic  idiosyncrasy  of  the  individual  analyst.  In  addition, 
the  scale  of  the  analysis  is  rather  variable:  where  the  data  are  dense, 
the  analysis  is  on  a  smaller  scale  than  where  sparser  data  necessitate  a 
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greater  degree  of  interpolation.  Also,  studies  baaed  upon  the  small 
sample  of  a  few  synoptic  situations  may  lack  generality  in  their  results. 

All  these  factors  combine  to  subject  synoptic  studies  of  vorticity  to  a 
large  indeterminate  degree  of  variance. 

*  It  is  proposed,  therefore,  that  in  this  study  objective  methods  of 
vorticity  computation  be  used  to  test  the  hypothesis  that  a  significant  dif¬ 
ference  exists  between  geostrophic  and  observed  vorticity  at  500 mb.  Such 
a  method,  though  sensitive  to  observational  error  and  unrepresentativeness, 
has  the  advantage  of  greater  sensitivity  to  real  differences.  Estimates  of 
observational  error  and  of  unrepresentativtuasss  can  be  made.  The  ulti¬ 
mate  question  is  whether  a  valid  estimate  of  real  differences  in  the  two 
vorticities  can  be  detected  above  this  background  level  of  error. 

STUDY  I 

1.  Methods 

In  Study  I  four  networks  of  five  rawinsonde  stations  each  were  used 
for  objective  calculations  at  500  mb  of  geostrophic  vorticity  from  the  re¬ 
ported  heights  and  of  observed  vorticity  from  the  components  of  the  re¬ 
ported  winds.  The  networks  are  shown  in  figure  1  and  the  individual 
stations  are  listed  in  Table  1. 

A  300  km  scale  was  initially  dosen  for  this  study,  on  the  basis  of 
its  rather  widespread  use.  Since  the  geostrophic  vorticity  is  computed 
from  the  second  derivative  of  the  height  field  and  the  observed  vorticity 


Fill.  M4I«hM«  MtMffct  iu<  far  ttrllilly  nantillHi 


10 


Table  1.  Networks  and  stations  used  for  computations. 


Network  I 

Network  II 

Network  III 

Network  IV 

SQ  (center) 

Washington 

Fort  Worth 

North  Platte 

Ogden 

(DCA) 

(FWH) 

(LBF) 

(HIF) 

Sj  ("west") 

Greensboro 

Big  Springs 

Denver 

Ely 

(GSO) 

(BGS) 

(LRY) 

(ELY) 

S2  ("north") 

Pittsburgh 

Oklahoma  City  Rapid  City 

Boise 

(PIT) 

(OKC) 

(RAP) 

(BOI) 

S3  ("east") 

Hempstead 

Shreveport 

Omaha 

Lander 

(HEM) 

(BAD) 

(OMA,  OFF) 

(LND) 

S4  ("south") 

Norfolk 

San  Antonio 

Dodge  City 

Grand  Junction 

(orfi; 

(SAT) 

(DDC) 

(GJT) 

from  the  first  derivative  of  the  wind  field,  the  grid  size  used  for  determin¬ 
ing  observed  vorticity  should  be  one-half  the  size  used  for  computing  geo- 
strophic  vorticity.  It  was  not  possible,  however,  to  locate  such  a  smaller 
network  of  wind  reporting  stations  within  the  larger  network  of  height 
reporting  stations.  It  was  assumed,  therefore,  that  the  magnitude  of  the 
wind  components  varied  linearly  between  the  central  station  and  each  of 
the  outer  stations,  an  asaumption  which  states  in  effect  that  observed 
vorticity  of  a  600  km  scale  is  equivalent  to  that  of  a  300  km  scale.  This 
assumption  of  linear  wind  change  is  commonly  used  to  interpolate  be¬ 
tween  data  on  synoptic  charts  and  is  the  best  estimate  of  wind  change 
that  can  usually  'be  made  {Rapp,  1952), 

Only  those  networks  were  chosen  for  this  study  whose  outside 
stations  were  approximately  300  km  from  the  central  station  and  so 
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located  that  perpendicular  axes  through  the  central  station  could  be 
drawn  close  to  each.  Only  the  grid  axes  of  network  II  could  be  oriented 
north-south;  the  axes  of  the  other  networks  were  rotated  so  as  to 
minimize  the  mean  distance  of  the  outside  stations  from  the  axes.  A  linear 
interpolation  process  was  then  used  to  transform  station  values  into  grid 
values  by  the  application  of  small  corrections.  As  an  example  of  how  this 
procedure  was  applied,  we  shall  only  consider  network  I,  since  the  other 
networks  were  treated  similarly. 

A  diagram  of  network  I  is  shown  in  figure  2.  The  axes  of  the  gird 
were  rotated  40  degrees  counterclockwise,  to  satisfy  the  condition  that 
the  outside  stations  be  the  shortest  mean  distance  from  the  axes.  Next, 
the  transformation  of  station  to  grid  values  was  developed  in  terms  of  a 
general  scalar  quantity  S  to  be  later  applied  specifically  to  the  heights 
and  the  wind  components.  To  obtain  grid  value  Sj' ,  for  example,  we 
first  interpolate  the  value  for  the  intersection  of  the  line  SjS2  with  the 
X'  axis.  The  interpolated  value,  Sj*,  can  be  expressed  as 

S,*Si  , 

Si*  =  sz  + 

^1^2, 

If  we  now  substitute  the  value  of  S^*S j/  Si  S2  , 

Sf  s  Sj  +  ,167(S2-  Sj) 
or 

SJ*  *  .833 S^  +  .167S2  . 
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The  value  of  S^'  can  now  be  obtained  by  interpolation  between  S*  and 
SQ  as  follows: 


Sj'  =  s*  -  (S_  -  Sj*)  . 

?!*S0 


If  we  now  substitute  the  value  of  ySj*S,'/S,*S  , 

Si  VS,*-  .  104(SQ  -  Sj*)  , 
or 

S,'  =  1. 104  S,*  -  .  104  SQ  . 

If  we  substitute  the  value  derived  for  Sj*  , 

S,  =  .92S,  +  .18S2  -  .10So  . 

In  a  similar  manner  the  values  of  82' »  S31  and  S4  can  be  obtained.  The 
procedure  can  be  summarized  as  follows: 

1.  Interpolate  linearly  to  obtain  the  value  at  the  intersection  of  an 
axis  ,  with  the  line  joining  two  stations  values. 

2.  When  the  intersection  is  outside  the  grid,  interpolate  between  the 
intersection  and  the  central  station  to  obtain  the  value  at  the 
grid  point. 

3.  When  the  intersection  is  within  the  grid,  extrapolate  the  gradient 
between  the  intersection  and  the  central  4tation  to  obtain  the 
value  at  the  grid  point. 

The  following  is  a  summary  of  the  relationships  thus  derived  for 


the  four  networks: 


Network  I 


Si'  =  .92 Sj.  +  .18S2  ’  .10So 
S21  =  .60  S2  4  .I8S3  4  .22  SQ 
s3'  =  .96S3  +  .09  S2  -  .05So 
s4'  =  I.00S4  -  .22  SL  +  .  22  SQ 

Network  II 

Sj '  =  .93  Sj  +  .16S2  -  -09So 
S2*  =  .83  S2  +  .O6S3  *  .11  SQ 
S3'  =  .99  S3  +  .06  S2  -  .05  SD 
S4'  =  _70S4  +  .23  S3  +  .07  SD 
Network  III 

=  1.11  Sx  +  .06 S2  -  .17S0 
S2'  =  .86  S2  +  .  14  S3 
S3 1  =  .  69  S3  +  .  20  S2  +  .11  S0 
S4'  =  .96 S4  +  .  14 S3  -  .10So 
Network  IV 

Sj'  =  .76Sj  +  .11S4  +  .  13  SQ 
S2'  =  .71  S2  +  .29  SQ 
S3'  =  .68S3  +  .  14  S2  +  .  18  SQ 
S4'  =  1.08  S4  +  .05  Sj  -  .  13  SQ 
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(3) 
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Since  L  =  300  x  103m 

gQ  =  9.81m  sec"^ 

and  1  ft  =  .30B m  , 

?  3  =  -^-(h^  +  h2  +  h3'  +  h4'  -  4h0)  x  10'9  sec-1  , 
where  h'  denotes  the  interpolated  geopotential  height  of  the  isobaric  sur¬ 
face  in  tens  of  feet. 

If  we  consider,  for  example,  network  I  where  the  latitude  of  the 
central  station  is  38.8°  N  and  f  =  .914  x  10"^  sec  *  . 

5g3i  =  .363  (hj  +  h2'  +  h3'  +  h4'  -  4hQ)  x  10‘5  sec-1  . 

If  we  express  the  grid  values  in  terms  of  the  station  values  of  height,  as 
derived  previously, 

hj  =  .92b!  +  .18h2  -  .10ho 
h2  =  .60h2  +  .18h3  +  .22  hQ 
h3'  =  .96h3  +  .09h2  -  .05ho 
h4  =  1.00h4  -  .22h1  +  .22h0 

Upon  substituting  these  values, 

£g3I  =  . 363(0. 70h1  +  0.87h2  +  1.14h3  +  1.00h4  -  3.71hQ)  x  10"5sec-1. 
Similarly  we  obtain  for  the  other  networks, 

^g3n  =  AZ1  <°-93hl  +  l-05h2  +  1.28h3  +  0.70h4  -  3.96h0)  x  lO'Ssec"1. 

*5g3III  =  *347(1. llhj  +  1.12h2  +  0.97h3  +  0.96h4  -  4.16hQ)  x  10"5  sec'1. 

-^3IV  *  •  346(0, 81hj  +  0.85h2  +  0.68h3  +  1.19h4  -  3.53h0)  x  10"® sec"1. 

Throughout  the  remainder  of  this  paper  all  vorticity  values  will  be  under¬ 

stood  to  be  in  units  of  10’®  sec'1. 
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2.  Data 

The  original  plan  of  Study  I  was  to  compute  geostrophic  and  observed 
vorticity  from  500 mb  data  for  the  entire  year  1953.  The  preliminary  re¬ 
sults  were  such,  however,  that  this  study  was  limited  to  0300Z  data  for  the 
period  January-June  1953.  The  data  were  obtained  from  the  Daily  Upper 
Air  Bulletins,  original  teletype  data,  and  from  microfilmed  copies  of 
WBAN  Form  33,  Summary  of  Constant  Pressure  Data  (obtained  from  the 
National  Weather  Records  Center,  USWB).  WBAN  33  would  have  been 
used  exclusively  in  this  study,  but  the  coding  of  wind  directions  to  16  points 
on  this  form  made  such  data  virtually  unusable  except  as  a  check  for  gross 
errors  in  the  transmitted  values  of  wind  direction,  which  are  coded  to  a 
36-point  wind  scale.  After  convers}diulo  appropriate  units,  however, 
such  data  were  occasionally  used  to  fill  in  gaps  where  the  data  were  other¬ 
wise  complete  for  a  network. 

The  height  values  from  WBAN  33,  which  are  listed  in  tens  of 
meters,  were  converted  to  units  of  ten  feet,  and  these  values  were  also 
used  as  a  check  on  the  transmitted  data  and  to  fill  in  occasional  gaps. 

This  checking  system  eliminated  many  errors  of  coding  and  transmission, 
so  that  the  data  used  in  the  objective  calculations  were  the  best  available. 

The  wind  data  were  those  reported  for  the  500 mb  surface  with  two 
exceptions:  for  Pittsburgh,  where  the  RAWINS  are  reported  separately 
from  the  RAOBS,  the  wind  was  interpolated  to  the  height  of  the  500  mb 
surface;  and  for  Omaha,  where  only  sparse  P IB  ADS  were  available. 
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RAWIN  data  from  Offutt  AFB  were  used  by  interpolating  the  wind  to  the 
RAOB  height  reported  by  Omaha. 

Since  the  objective  system  of  calculation  required  five  wind  and 
height  reports  for  each  network,  the  rate  of  attrition  was.  quite  large. 
Table  2  lists  the  monthly  totals  for  each  network  where  the  data  were 
complete. 


Table 

2. 

Number  of  cases  computed. 

Network 

Jan. 

Feb.  Mar. 

Apr. 

May 

June 

Total 

I 

8 

8 

17 

12 

15 

17 

77 

II 

6 

12 

10 

14 

21 

25 

88 

in 

5 

2 

6 

10 

13 

19 

55 

IV 

6 

6 

12 

8 

3 

11 

46 

Total 

25 

28 

45 

44 

.Si.  . 

72 

266 

Statisti  cs 


Scatter  diagrams  (figs.  3,  4,  5,  6)  represent  graphically  the  relations 

between  ,  and  which  were  found  for  each  network.  The  note- 

ob  g  j 

worthy  features  of  these  diagrams  are  the  considerably  greater  range  of 

values  of  o  i  and  the  weak  relation  between  the  two  vorticities.  Of 

g3 

special  interest  are  the  large  negative  values  of  geostrophic  vorticity 
found  in  networks  I  and  II,  many  of  which  exceed  the  value  of  the  Coriolis 
parameter. 


The  following  statistics  were  computed;  the  arithmetic  means, 


FIG. 3.  Scatter  diagram  of  £ 


FIG-6-  Scatter  diagram  ot  £06  v#  Cg3  *0f  nt^ork3Z  . 


24 


5*3  and  ;  the  standard  deviations,  cr,^  and  or..  ^  ;  and  the  li_. ear 

correlation  coefficients,  r,  between  %  ->  and  S  ,  for  the  networks 

g3  06 

individually  and  collectively.  These  values  are  presented  in  Table  3 
below. 

Table  3.  Statistical  summary  of  the  comparative 

analysis  of  300km  geostr.ophic  vorticity 

( <Sa-i )  and  600km  observed  vorticity  ($”,). 
a J  06 


Network 

n 

5k3 

§06 

V 

°o6 

r 

I 

77 

-5.13 

-0.44 

7.18 

3.03 

.540 

II 

83 

-1.99 

-0.82 

10.22 

2.14 

.246 

III 

55 

7.89 

-0.73 

9.02 

3.37 

..323 

IY 

46 

1.52 

-1.37 

6.30 

1.93 

.211 

Total 

266 

-0.25 

-0.79 

9.76 

2.69 

.270 

4.  Conclusions 

The  linear  coefficients  are  so  low  as  to  be  insignificant.  The 
standard  deviation  of  ,  is,  for  the  combined  cases,  more  than 


g3 


three  times  larger  than  that  of  .5^.  This  striking  difference  in  the 
distributions  of  the  two  vorticities  can  be  attributed  to  a  combination 
of  the  following  factors:  observational  errors,  to  which  the  geostrophic 
vorticity,  computed  from  the  second  derivative  of  the  height  field,  is 
especially  sensitive  due  to  the  lack  of  smoothing:  unrepresentativeness, 
to  which  the  observed  vorticity  is  especially  vulnerable,  and  which  can 
be  reasonably  expected  to  flatten  out  the  range  of  values;  the  scale 
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effect,  which  was  introduced  into  this  study  by  the  assumption  that  the  wind 
components  varied  linearly  from  the  central  station  to  the  outkr. .  stations; 
and  finally,  any  real  differences  between  the  two  vorticities. 

Of  the  factors  mentioned  above  as  accounting  for  the  large  difference 

in  the  standard  deviations  of  5  _  and  estimates  of  the  effects  of 

g3  06 

observational  errors  and  unrepresentativeness  can  be  made  (and  will  be 
presented  later  in  this  paper).  The  real  differences  between  the  two  vorti- 
city  distributions  are  most  probably  obscured  by  the  systematic  scale 
effect  introduced  into  this  study.  If  the  assumption  of  linear  wind  varia¬ 
tions  was  not  valid,  the  original  hypothesis  was  not  subjected  to  a  fair 
test.  The  results  of  this  study  raise,  therefore,  an  important  question  as 
to  the  limitation  imposed  by  the  density  of  data  upon  the  scale  chosen  for 
vorticity  computations.  Interpolation  too  far  within  the  data  network  may 
seriously  affect  vorticity  computations  purporting  to  be  finer  than  the 
coarseness  of  the  data  will  allow.  Such  may  be  the  case  here,  where  it 
is  assumed  that  the  observed  vorticity  at  a  600km  scale  may  be  con¬ 
sidered  equivalent  to  that  on  a  300km  scale. 

Landers  (1955)  has  presented  an  example  of  the  effect  produced  on 
vorticity  computations  by  the  scale  of  wind  observations.  From  a  large- 
scale  analysis  of  synopitcally  reported  winds  at  30,000  feet,  anticyclonic 
relative  vorticity  was  computed  over  a  large  region, whereas  from  a 
small-scale  analysis  of  B-29  winds  at  the  same  level  for  approximately 
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the  same  time,  cyclonic  relative  vorticity  was  predominant  over  the 
same  region.  The  ratio  of  small  to  large  scale  maximum  gradients 
of  vorticity  was  six  to  one.  It  would  thus  appear  that  interpolation 
of  synoptic  wind  data  for  computing  smaller-scale  vorticity  may  be  sub¬ 
ject  to  large  errors. 
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STUDY  II 

l.  Methods 

The  results  of  the  previous  study  indicated  the  probable  seriousness 
of  the  assumption  of  linear  wind  variation  between  observing  stations.  It 
was  decided,  therefore,  to  alter  the  plan  of  study  so  that  the  effect  of  scale 
on  vorticity  computations  could  both  be  removed  and  investigated  separately. 
A  method  was  designed  to  compare  600  km  geostrophic  and  observed  vorti¬ 
city,  and  in  addition,  to  compare  300  km  and  600  km  geostrophic  vorticity. 

A  network  for  computing  600km  geostrophic  vorticity  was  chosen  whose 
central  station  coincided  with  that  of  network  III.  From  these  concentric 
networks  three  sets  of  vorticity  values  could  be  computed. 

The  larger  network  for  computing  600  km  geostrophic  vorticity, 
which  we  shall  designate  at  network  III',  is  shown  in  figure  7.  The  indi¬ 
vidual  stations  are  listed  in  Table  4.  The  axes  of  network  HI',  chosen 


Table  4.  Stations  used  for  computing  600  km  geostrophic  vorticity. 


Station 
SQ  (centerl 
Sj  ("west") 

52  ("north") 

53  ("east") 

54  ("south") 


Network  III' 
North  Platte  (LBF) 
Lander  (LND) 
Bismarck  (BIS) 
Columbia  (CBI) 
Amarillo  (AMA) 
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to  minimize  the  interpolation  necessary  for  computing  grid  values  of  height 

from  station  values  was  plareed  clockwise  from  the  axes  of  network  III. 

The  relationships  previously  derived  in  Study  I  for  ^  and 

"C  were  again  applied  to  this  study;  the  same  methods  used  in  their 

g3HI 

derivation  were  applied  to  the  larger  network  and  the  following  expression 
thus  obtained  for  computing  the  600  km  geostrophic  vorticity: 

S  g6mi  =  0.087  (0.86h1  +  0.80h2  +  1.08h3  +  1.13h4  -  3.87hQ)  x  10‘5sec_1. 
Three  vorticity  values  were  then  computed:  ^g3  »  0b  *  anc*  ^g6 
subscript  denoting  the  network  is  no  longer  necessary)*  Only  those  cases 
were  selected  where  all  three  values  could  be  obtained,  so  that  the  com¬ 
parison  between  and  ^  and  that  between  ^3  and  were  based 

upon  the  same  observations. 

2.  Data 


A  total  of  385  cases  were  selected  from  0300Z  and  1500Z  upper  air 
data  for  the  year  1953.  The  following  is  the  monthly  distribution  of  cases 
for  the  months  January  to  December,  respectively:  18,  14,  15,  36,  27, 

39,  47,  47,  44,  40,  31,  27.  It  is  apparent  that  the  summer  months  have  a 
larger  number  of  cases  than  the  winter  months.  In  an  attempt  to  eliminate 
some  of  the  bias  resulting  from  the  unequal  seasonal  distribution,  the  cases 
were  divided  into  two  groups;  the  "winter"  group  of  months  were  Jan. , 

Feb. ,  Mar.  ,  April,  Nov.  ,  and  Dec. ,  and  the  "summer"  group  of  the  months 
were  May,  June,  July,  Aug.,  Sept.,  and  Oct.  The  winter  cases  total  141; 


the  summer  cases  244. 
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3.  Statistics 

V 

A.  The  relationship  between  $  ^  and  $0£* 

The  relationship  between  and  >^o6  *s  represented  separately 

for  winter  and  summer  cases  by  scatter  diagrams  (figs.  8  and  9).  Table  5 
summarizes  the  following  statistics:  the  number  of  cases,  n,  the  arith¬ 
metic  means,  anc*  ~'o6*  t^e  standard  deviations,  and  <tq£  ;  the 

linear  correlation  coefficients,  r  ;  the  standard  error  of  estimating  ^  / 
from  5 ;  the  standard  error  of  estimating  JTq£  from-^g^,  Sq^  ; 
and  the  root  mean  square  of  the  absolute  differences  between  Sq£)  , 


^g6’ 


rms 


3V 


Table  5.  Statistical  summary  of  comparative  analysis 

of  600  km  geostrophic  and  observed  vortitities. 


Group 

n 

5g6 

^6 

V 

°o6 

r 

S«6 

So6 

rmsj6 

Winter 

141 

0.97 

-0.33 

3.55 

2.98 

.764 

2.28 

1.92 

2.65 

Summer 

244 

1.16 

-1.25 

2.69 

2.60 

.653 

2.04 

1.97 

.  3.  26 

Total 

385 

1.09 

-0.91 

3.03 

2.78 

.689 

2.20 

2.02 

3.05 

There  is  some  evidence  of  a  seasonal  effect  in  the  fact  that  for  the 
summer  group  the  correlation  coefficient  is  somewhat  lower,  the  spread 
of  the  means  is  slightly  increased  and  the  standard  deviations  of  both 
vorticities  are  decreased.  The  standard  error  of  estimating  5  ^  from 
^  is  slightly  lower  for  the  summer  group;  the  root  mean  square  of  the 
absolute  differences  is  larger  for  the  summer  cases.  If  we  consider  all 
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the  cases,  the  separation  of  the  means  of  the  vorticities  is  equal  to  slightly 
over  20  percent  of  the  value  of  the  Coriolis  parameter;  the  standard  errors 
of  estimating  one  vorticity  from  the  other  is  also  about  20  percent  of  the 
Coriolis  parameter;  the  root  mean  square  of  the  absolute  differences  is  over 
30  percent  of  the  value  of  the  Coriolis  parameter.  The  correlation  co¬ 
efficient  for  all  the  cases  indicates  that  47  percent  of  the  variance  of  one 
vorticity  can  be  accounted  for  by  the  other,  —for  the  winter  cases  alone,  58 
percent;  for  the  summer  cases  alone,  43  percent.  The  following  is  a  list 
of  the  equations  of  the  regression  lines: 

1.  All  cases  :  S  =  0.75  +  1.77 

^'ob  -0.63^-1.60 

2.  Winter  cases  :  g£>  =  0.91  +  1.27 

$'o(>  =  0.64  5g6- 0.95 

3.  Summer  cases:  gS  =  0.67-^k  +  2.00 

J’o6  =  0.63  Sg6-  1.98 

If  we  assume  typical  values  of  ^ and  use  the  regression  equations 

given  above  to  estimatedvKhlues  of  %  ,  we  obtain  the  values  shown  in 

go 


Table  6. 
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Table  6.  600  km  geostrophic  vorticities  estimated 

from  regression  equation  for  various 
values  of  600km  observed  vorticity. 


winter 


■10 

-  7.8 

5 

-  3.3 

0 

1.3 

5 

5.8 

10 

10.4 

15 

14.  9 

summer  all  cases 


4.7 

-  5.7 

1.4 

-  2.0 

2.0 

1.8 

5.4 

5.5 

8.  7 

9.3 

12.1 

13.0 

If  we  assume  typical  values  of  and  similarly  compute  the 

estimated  values  of  J  ,,  we  obtain  the  values  shown  in  Table  7. 

oo 


Table  7.  600  km  "observed"  vorticities  estimated 

from  regression  equation  for  various 
values  of  600  km  geostrophic  vorticity. 


-10  -  7.4 

-5  -  4.2 

0  -  1.0 

5  2.2 

10  5.4 


summer 

all  cases 

-  8.3 

-  7.9 

-  5.1 

-  4.8 

-  2.0 

-  1.6 

1.2 

1.6 

4.3 

4.7 

B.  The  relationship  between  and  % 

g6  g3 

The  relation  between  an<^  ^3’  geostrophic  vorti- 

cities  on  the  two  different  scales,  is  represented  for  the  separate 

groups  of  winter  and  summer  cases  by  scatter  diagrams  (figs.  10  and 

11).  Table  8  represents  statistics  similar  to  those  listed  in  Table  5 

for  and  where  now  y>g£  is  the  standard  error  of  estimating 

from  ~^g3»  and  Sgj  is  the  standard  error  of  estimating 

from  o  l  • 

g& 


Table  8.  Statistical  summary  of  comparative  analysis  of 
300  km  and  600  km  geostrophic  vorticities. 


Group 

n 

l 

V 

r 

s' , 
g6. 

JjLl 

Winter 

141 

0.92 

5.52 

3,55 

7.16 

.726 

2.43 

4.93 

Summer 

244 

1.16 

7.44 

469 

7.44 

.740 

1.81 

5.01 

Total 

385 

1.09 

6.74 

3.03 

7.39 

.723 

2.09 

5.11 

The  evidence  of  a  seasonal  effect  is  not  as  strong  as  was  found 
in  the  study  of  ^g^  and  The  correlation  coefficients  for  the 

separate  winter  and  summer  cases  show  remarkable  agreement;  how¬ 
ever,  the  spread  of  the  means,  as  before,  is  greater  for  the  summer 
group.  The  winter  and  summer  values  of  and  are  in  close 
agreement.  The  standard  error  of  estimating  5?  ^  from  is 

very  large,  approximately  55  percent  of  the  value  of  the  Coriolis 


parameter. 


FI6.II.  Scattor  diagram  of  C-6vs.  £a,  t  Summtr  coots. 
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If  we  consider  all  the  cases,  the  differences  in  the  means  and  the 
standard  deviations  are  quite  striking.  Yet,  the  two  vorticities  are 
rather  highly  correlated,  since  52  percent  of  the  variance  of  one  can 
be  accounted  for  by  the  other.  The  following  is  a  list  of  the  equations  of 
the  regression  lines: 


All  cases  : 

S'. 

g*> 

0.30 

Sg3 

-  0.91 

_ 

°  g3  ' 

1.76 

Sg6 

+  4.82 

Winter  cases  : 

0.36 

-^g3 

«  1.02 

^g3‘ 

1.47 

+  4.09 

Summer  cases: 

0.27 

V 

-  0.82 

2.05 

Sg6 

+  5.06 

Le  typical  values 

of  Sg3  ‘ 

and  use  the  regression 

as  given  above  to  compute  the  estimated  values  of  we  obtain  the 

values  shown  in  Table  9. 

If  we  assume  typical  values  of  arid  compute  the  estimated 

values  of  5  ,  ,  we  obtain  the  values  in  Table  10. 
g3 

If  one  were  to  assume  a  linear  variation  of  the  height  differences, 

(hj  -  hQ),  etc.,  from  one  grid  scale  to  another,  one  would  expect  a  given 
value  of  to  be  one-half  that  of  the  corresponding  value  of  ,  since 

the  computation  of  geostrophic  vorticity  depends  inversely  upon  L^,  the 
square  of  the  grid  differentiation  difctaftce.  It  would  appear,  from  the 
sample  calculations  made  from  the  regression  equations,  that  the  use 
of  positive  values  of  to  calculate  produces  an  estimate  of  the 
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Table  9.  600km  geostrophic  vorticities  estimated 

from  regression  equation  for  various 
values  of  300  km  geostrophic  vorticity. 


Table  10.  300  km  geostrophic  vorticities  estimated 

from  regression  equation  for  various 
values  of  600  km  geostrophic  vorticity. 
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larger-scale  geostrophic  vorticity  considerably  less  than  one-half  the 
value  of  the  smaller-scale  vorticity.  On  the  other  hand,  the  use  of 

to  estimate  ^ produces  values  of  the  smaller-scale  vorticity 
which  are,  generally  speaking,  about  (twice  that  of  thq  larger  scale 
vorticity. 

4.  Discussion  of  Errors  and  Unrepresentativeness 

Before  attempting  a  final  evaluation  of  the  statistical  results,  an 
estimate  should  first  be  made  of  the  variance  introduced  into  the  calcu¬ 
lations  by  errors  in  the  measurement  of  heights  and  winds  and  by  small- 
scale  fluctuations.  We  shall  consider  here  the  error  standard  deviation. 
From  equation  (3),  it  is  obvious  that  the  calculation  of  geostrophic  vorti¬ 
city  is  most  sensitive  to  errors  in  measuring  the  height  of  the  central 


station.  The  error  standard  deviation  of  the  geostrophic  vorticities  can 
be  shown  (Scarborough,  1930)  to  be 


If  we  assume  that  the  error  variances  of  the  height  values  are  equal  for 
the  five  stations,  we  obtain,  after  simplifying 

fZo  Ofcg 

- 2 -  (5) 

fir 

If  we  now  use  Rapp's  (1952)  estimate  of  the  error  variance  of  the 

2  2 

height  of  the  500 mb  surface,  =  104.5m  ,  we  obtain  as  the  error 
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standard  deviation  of:  ^  _  , 

°g3 

=  5.2  x  10"^  sec  *  . 

Similarly,  we  obtain  as  the  error  standard  deviation  of  ^  ,  , 

go 

a  value  one-fourth  that  of  /cr  „ 

3g3 

If  we  now  also  consider  equation  (1)  used  in  computing  the  observed 
vorticity,  we  can  express  the  error  standard  deviation  of  the  observed 
vorticity  as 


a  JAX.zJA 

$°S/  \  9v,  )  v3  1  3v, 


\8vl)  vl  \8u2 


*£)  "u2  +  \ K  )  ru4 


2  2  2  2 

If  we  assume  that  <r  ,  *  e  .  and  <r  =  <r  .  ,  it  can  be  shown  that 

v3  vl  u2  u4 


From  Rapp's  results  the  average  error  variances  of  the  east-west 
and  north-south  wind  components  between  5.5  and  6.0km  above  sea 

level  (the  approximate  range  of  the  height  of  the  500  mb  surface)  are 

Z  2  -2  Z  Z  .7 

found  to  be  <r  =  .65m  sec  and  <r  =  .79m  sec*  .  Thus  the  error 
u  v 

standard  deviation  of  S  .  is  found  to  be 

ob 

~  *24  x  10"^ sec"*  . 

These  estimates  of  the  error  standard  deviations  must  be  viewed,  however, 
as  being  on  the  low  side,  since  Rapp's  statistical  values  were  based  upon 
data  more  accurate  than  that  obtained  in  daily  observational  routine. 


In  addition  to  the  variance  produced  by  errors  in  measuring 
heights  and  winds,  the  variance  of  small-scale  fluctuations  should  be 
considered.  It  is  probable  that  this  variance  is  of  considerable  im- 
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portance,  since  the  wind  observations  are,  in  effect,  non-simultaneous,  and 
the  small-scale  fluctuations  are  of  a  short  enough  period  (5  minutes  to 
3  hours}  to  seriously  affect  vorticity  computations  based  on  the  reported 
winds  of  five  stations. 

Rapp  lists  as  possibly  useful  approximations  to  the  average  magni- 

A  2 

tudes  of  the  variances  of  the  small-scale  fluctuations,  values  of  o-^,  and 
rs2 

,  for  various  heights  as  determined  from  four  separate  experiments. 
The  values  which  he  obtained  for  o-^,  and  o’  ,  between  5  and  6  km  were 
then  averaged.  The  mean  of  seven  values  of^^,  was  3.3m^sec'^;  the 
mean  of  five  values  of  was  lO.lm^sec"^  (due  to  the  occurrence  of 
two  values  exceeding  20m^sec"^).  If  we  assume  that^^,  and  o\^,  are 
each  equal  to  a  value  of  approximately  5m^sec"^,  the  standard  devi¬ 
ation  of  the  observed  vorticity  produced  by  small-scale  wind  fluctuations 
is 

<Tgo6  s  1*7  x  10“^ sec” *  . 

The  above  values  of  the  standard  deviations  of  vorticity  resulting  from 
observational  error  and  small-scale  fluctuations  are  proposed  as  rough 
qualitative  estimates.  It  would  seem  reasonable,  however,  to  assume 
that  the  variance  of  observed  vorticity  due  to  small-scale  wind  fluctu¬ 
ations  significantly  exceeds  that  due  to  errors  in  wind  measurement. 


and  that  their  combined  variance  is  comparable  to  that  of  the  geostrophic 
vorticity  of  the  same  scale. 

5.  Conclusions 

On  the  basis  of  the  results  of  this  study,  it  is  concluded  that  geo¬ 
strophic  and  observed  relative  vorticity,  when  measured  on  t’  e  same 
scale  of  600  km,  are  closely  related,  especially  during  the  winter.  In 
the  "winter"  months  58  percent  of  the  variance  of  one  vorticity  can  oe 
accounted  for  by  the  other;  in  the  "summer"  months  only  43  percent  can 
be  accounted  for.  No  great  disparity  could  be  found  in  the  standard  de¬ 
viations  or  in  the  standard  error  of  estimating  one  vorticity  from  the 
other,  except  for  a  slight  tendency  for  to  be  larger  than  o^o£, 

an  effect  more  pronounced  in  the  winter  cases.  The  geostrophic  and 
observed  vorticities,  although  closely  related,  are  not  identical  in 
magnitude.  The  root  mean  square  values  of  the  absolute  differences 

06 

nitudes,  even  if  we  consider  the  effect  of  observational  errors  and  un¬ 
representativeness.  The  standard  error  of  estimating  one  vorticity 
from  the  other  is  considerable,  about  20  percent  of  the  value  of  the 
Coriolis  parameter  (for  the  region  of  study,  f  ^  9  x  10  sec"  h  If  we 
consider,  moreover,  the  variance  of  small-scale  wind  deviations  in  ad¬ 
dition  to  the  error  variance  of  wind  measurement,  it  is  probable  that 
observed  vorticity  cannot  be  more  accurately  measured  than  the 


indicate  rather  significant  differences  in  their  mag- 


geostrophic  vorticity. 
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The  results  of  the  study  of  300  km  and  600  km  geostrophic  vor.tii 
city  indicate  the  important  effect  of  grid  size  upon  geostrophic  compu¬ 
tations.  This  effect  is  mainly  brought  about  by  second  order  differen¬ 
tiation  of  the  height  field,  so  that  even  with  a  linear  variation  of  height 
differences,  doubling  the  grid  size  reduces  the  relative  vorticity  by  one- 
half.  The  results  seem  to  show  that  the  vorticity  centers  are  of  such 
small  scale  that  doubling  the  grid  size  from  300km  to  600km  will  re¬ 
duce  the  values  of  vorticity  even  further.  The  effect  produced  by  errors 
in  height  measurement  upon  the  computation  of  geostrophic  vorticity  is 
also  dependent  on  scale  since  the  effect  of  doubling  the  grid  size  is  to 
reduce  error  variance  to  one-fourth  that  of  the  smaller  scale. 

The  first  study  of  and  failed  because  it  was  based  upon 

a  fallacious  assumption  that  vorticities  of  different  scales  were  com¬ 
parable.  The  scale  of  measurement  has  such  an  important  effect  upon 
the  computation  of  geostrophic  vorticity,  that  it  might  also  be  expected 
to  have  a  similar  effect  upon  observed  vorticity,  although  this  effect 
could  not  be  directly  evaluated  in  this  study. 
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